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By E. R. G. Eckert^ John N. B. Livingood^ and Ernst I. Prasse 


SUMMARY 

In transpiration cooling of various structural elements in gas tur- 
bines^ the coolant has to be ducted within passages to the porous walls 
through which it is ejected into the gas stream. The passages, often 
arranged in rotating parts, have to be designed in such a way as to en- 
sure the proper local distribution of the coolant. In this report, a 
method is presented by which either the local permeability necessary for 
a prescribed distribution of the coolant flow or the coolant-flow dis- 
tribution resulting from a prescribed local permeability can be predicted. 
The method is based on a one- dimensional treatment of the gas flow through 
a rotating channel with varying cross section and partially porous walls. 
The inlet pressure into the channel and the outside pressure along it are 
assumed prescribed. It is also stipulated that the passage ends blindly. 
However, the method can easily be extended to cover the situation where 
a certain mass flow leaves the open end of the passage. 

The method was applied to a determination of conditions in rotating 
turbine blades with transpiration- cooled walls. The cooling air is as- 
sumed to flow from the blade root through several channels to the porous 
skin of the blades. For prescribed coolant-flow ejection rates necessary 
to maintain a constant porous wall temperature, considerable reduction 
in wall permeability from passage entrance to passage tip (blade root to 
blade tip) is required. In fact, such required variations will be ex- 
tremely difficult to fabricate, and compromises between prescribed wall 
temperature and coolant-flow ejection may be necessary. For prescribed 
locally constant permeability, the mass-flow ejection rate increased from 
blade root to blade tip. The relative increase along the blade became 
smaller when the inlet pressure of the coolant at the blade root was in- 
creased. A large inlet pressure is therefore conducive to a uniform flow 
ejection rate and, accordingly, a more uniform blade wall temperature. 

In the investigated range, a passage area variation had practically no 
effect on flow ejection rates. For low Mach numbers at the passage inlet 
(below approximately 0.2 in the examples), the internal pressure distri- 
bution may well be approximated by consideration of rotational effects 
only. This leads to a considerable simplification in the calculation 
procedure . 


2 


NACA TN 3408 


INTRODUCTION 

Transpiration cooling is considered as an effective means of keeping 
different structural elements in gas turbines at a temperature level suf- 
ficiently low for the material from which the elements are fabricated. 
Application of this cooling method is considered^ for instance^ for cool- 
ing the walls of combustion chambers and of rotating turbine blades. In 
both cases ^ passages have to be provided within which the coolant is 
ducted to the porous surfaces before it is finally ejected into the gas 
stream. 

When transpiration cooling is applied to the rotating blades of the 
gas turbine^ these passages will have to be arranged in the interior of 
the blade, for instance, as indicated in figure 1. The passages are 
closed at the top of the blade, and the cooling air leaves the passages 
through the porous blade surfaces. The flow of the cooling air through 
each of the passages is subject to friction forces along the passage 
wall, to centrifugal forces that tend to increase the pressure from the 
blade root towards the blade tip, and to pressure changes connected with 
the change in momentum of the cooling air on its way aJLong the passage 
interior. The problem arises as to what local permeabilities along the 
passage are required to ensure a certain flow distribution or what the 
local flow distribution of the cooling air will be as it passes the po- 
rous surface with prescribed permeability and is ejected into the gas 
stream. In both cases, the calculation is complicated by the fact that 
the pressure outside the blade varies locally in radial direction as well 
as around the blade circumference. 

A proper design of the blade requires the ability to predict the 
local flow rates through the porous blade walls. It is the purpose of 
this report to present a method by which the flow through passages with 
porous walls of the form indicated in figure 1 may be calculated. The 
method is based on a one- dimensional treatment of the flow and is essen- 
tially the same as presented in references 1 and 2. It is, however, 
especially adapted to the problem under investigation. The basic equa- 
tions necessary for a solution of the problem are derived at first in a 
general way, so that they can be applied generally to calculate flow 
conditions in rotating passages. These equations are then adapted to 
the speciaJ- conditions expected in gas turbine blades. The calculation 
procedure is discussed in detail, and the importance of the various pa- 
rameters affecting the flow distribution is investigated in a niomber of 
numerical solutions. 

SYMBOLS 

The following symbols, with consistent units, are used: 

A cross-sectional (flow) area of passage 
a slope of dimensionless temperature ratio, t = 1 + a§ 
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b 

length of porous section of passage measured in circumf erentieuL 
direction 


quantity containing permeability coefficient (see eq. (ll)) 

^K,dim 

dimensionless quantity containing permeability coefficient 
(see eq. (l7)) 

w ^1 

cn 

viscous-resistance coefficient 

C\D 

r 

^2 

inertial-resistance coefficient 

"P 

specific heat 

Dh 

hydraulic diameter of passage (= circumference) 

o f 

0) 

friction coefficient 

o 

V ^ 

function 

'Si 

^ K 

permeability coefficient (see eq. (9)) 

' L 

length of passage 

M 

Mach number at passage inlet 

m 

mass coolant flow 

N 

CDL/Ay-yRTr 

W” 

^Pi,r/|RTr (Va)^ 

n 

exponent 

P 

pressure 

R 

gas constant 

r 

radius to element measured in plane of rotation (see fig. l) 

T 

temperature of coolant 

T 

• g 

effective gas temperature 

t 

thickness of porous wall 
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V flow velocity throiigh porous wall (based on total porous surface 

area) 

W velocity of main flow through passage 

X distance from passage entrance (see fig. l) 

a dimensionless area ratio^ 

p angle between velocity vector and direction of increasing radius 
X ratio of specific heats ^ 1.4 

|jL dimensionless coolant mass velocity ratio^ (m/A)/(m/A)p 

V kinematic viscosity 

4 dimensionless distance from passage entrance^ x/L 

iTg dimensionless external pressiire ratio, Pe/Pi^r 

dimensionless internal pressure ratio, Pi/Pi^r 

p dens ity 

T dimensionless temperature ratio, T^Tp, except for Ty (shearing 
stress at wall) 

03 angular velocity 

Subscripts: 

e external 

i internal 

r root (passage entrance) 

s considers rotating effects only 

w wall 

0 NACA standard conditions 

Superscript: 

* total condition 
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AMLYSIS 


Derivation of Basic Equations 

The basic equations used in this report are developed with the help 
of figure 1. This figure shows a portion of a channel^ one wall of which 
is porous. The channel may rotate around the axis 1-1 with the angular 
velocity o). Attention is fixed on a portion of fluid between the cross 
sections A and A + dA. This region is located at the distance r 
from the axis 1-1^ and the axis of the channel is inclined under an angle 
P against the radial direction^ where p is measured between the ve- 
locity vector and the direction of increasing r. The flow through the 
passage is assumed to be steady relative to the channel. 


For steady-flow problems^ application of Newton* s second law of mo- 
tion may be made by equating the net force acting on a control surface 
and the body forces acting on the fluid particles within the control area 
to the increase in momentum of the stream flowing through the stationary 
control surface. In the present case this law must be applied to a con- 
trol area that is at rest relative to the channel. In figure 1 such a 
control area is indicated bounded by the two cross sections A and 
A + dA and by a surface located in immediate proximity to the inside 
wall of the channel between the two cross sections. Forces acting in 
flow direction on the control surface will be pressure forces in the two 
cross sections and along the channel wall^ friction forces acting eilong 
the channel wall in the portion between the two cross sections^ centrifu- 
gal forces^ and Coriolis forces. The sum of all pressure forces is 


P±A - (Pi + iiPi)CA + dA.) + 



dA 


if the pressure on the surface bounding the channel walls is approximated 
by pj^ + dpj^/2. Neglecting second-order terms reduces the pressure 

forces to 


- A dpjL 


When the friction factor f is defined by the equation 


T 


w 


8 


(t.^^ denoting the sheaxing stress at the wall), the friction force can 
be expressed in the following way: 
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Body forces^ caused by centrifugal acceleration of the fluid particles^ 
can be written as 


2 

pAro) cos p dx 

The Coriolis force is oriented normal to the flow direction through the 
passage and has no component in the direction of the passage axis. The 
increase in momentum of the stream flowing through the control surfaces 
is 


or 


(m + dm)(W + dW) - mW 


in 


m dW + W dm 

Equating the various force terms with the change in momentum results 


A dpi - ~ ^ pAro)^ cos pdx = mdW + Wdm (l) 

2 


Use of the continuity equation 


m = pAW 


the equation of state 


Pi = pRT 


( 2 ) 

( 3 ) 


and the relation 


dm 


it) 


m dA 


dx dx ‘ A dx 


in equation (l) leads to 


‘ ■ ©■ 3 


(!) 


Pi 2 a 2RT /m 

= (rr: rcu cos p - ( — i 

dx RT Pi \Ay ix 


m\- I R \/dT ^ T ^ ^ 


a) ypWldx A dx 2D^j 


( 4 ) 


( 5 ) 
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The static temperature T is determined by the rate at which internal 
energy in the flowing gas stream is converted into kinetic energy^ as 
well as by the rate with which heat is added to or subtracted from the 
gas stream. For high-velocity flow in which the conversion of internal 
energy into kinetic energy becomes an important factor^ it woiild be more 
appropriate to calcijlate with the total temperature^ which is directly 
connected with the rate of heat addition^ than with the static tempera- 
ture. This conversion can be made with the help of the following 
equation: 


The change from the static to the total temperature^ however^ complicates 
equation (S) considerably. On the one hand^ the flow velocities in the 
passages are usually small enough to make the temperature change connected 
with the conversion of internal into kinetic energy small. On the other 
hand^ the rate of heat addition to the coolant on its way through the 
flow channel usually is not too well known and has to be estimated. For 
these reasons^ the static temperatirre is maintained in the following 
calculations. It is assumed that the vaxiation of the static temperature 
along the coolant flow channel is known. 

In equation (5)^ the specific mass flow m/A^ the internal' pressure 
p^^ the radius r^ the tangential velocity ra>^ and the cross-sectional 

area A are generally functions of the distance x from the channel 
entrance measured along the channel axis. The temperature T may also 
vary along x. 

A second expression including the same variables as functions of 
the distance from the passage entrance can be obtained from a considera- 
tion of the pressure drop throiigh the porous wall when the pressure dis- 
tribution along the outside of the passage is prescribed. The mass 
velocity of the coolant ejected through the porous wall at any location 
multiplied by the width of the porous portion of the channel wall meas- 
ured in circumferential direction equads the rate of change of the cool- 
ant mass flow passing through the passage at that locationj that is^ 


The quantity pV is the average mass-flow rate over the channel width 
b. The negative sign appears in the left member of equation (6) when 
the velocity is considered positive for flow leaving the channel throiigh 
the porous wall. Reference 3 shows that the mass velocity of a gas 
flowing through a porous wall is connected with the difference in the 
squares of the pressures acting on both sides of the wall^ through the 
equation 




( 6 ) 


8 


NACA TN 3408 


2 2 

Pj " P o 

= Ci(2ETpv)(pV) + Cg 2RT(pV)^ (?) 

Values of and must be determined experimentally for each porous 

material. 


Eqxiation (7) may b6 referred to measurements at standard conditions 
by transformation into (see ref. 4) 


P 


2 

i 


t 



e 




( 8 ) 


The parameters appearing on both sides of equation (8)^ when plotted on 
log-log coordinates for several types of porous materials, indicate an 
almost linear variation over a fairly wide range of flow for any of the 
materials. As a consequence, equation (8) may be well approximated by 


pV 


pq^o 

pv 


K / PQ'^O X^ ^ 
t pv y T 




( 9 ) 


where K is the permeability of the material. On a log-log plot of 


p? - p^ 

■*^1 -^e 


'PoV ^ 


pv 


against pV(pqVq/pv), the intercept of the curve for 
a given porous material on the log pV(pqVq/pv) ajcis is the value of n 


log K for that material. This definition of K is similar, but not 
identical, to that normally used. For the purpose of this calculation, 
equation (9) is simplified to 


pV = Ck(p| - 4 )“ 


( 10 ) 


by introduction of the parameter 


Ck - 



( 11 ) 


The values of Crr and n are functions of the porous material^ they 
also vary slightly with the specific mass flow pV(poVo/pv). The value 
of is a function of the material permeability and thickness, and of 

cooling-air properties based on porous-wall temperature. Substitution 
of equations (4) and (lO) into eq-uatlon (6) gives the desired relation: 
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dx 


m ^ 

A dx 




( 12 ) 


Equations (5) and (12) are two equations from which two unknown values 
can be calculated. 


In order to make the results of calculations more generally appli- 
cable^ it is advantageous to change equations (5) and (12) to dimension- 
less form. For this purpose^ the following substitutions are made (the 
dimensionless variables are (^enoted by Greek letters corresponding to 
the English letters for the dimensioned quantities^ and the constants by 
English capital letters): 


Reduction to Dimensionless Form 




Pi 



^ Pi,r 



(13) 






J 
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It can be shown easily that the parameter M is the Mach number of the 
flow at the passage inlet. The parameter N" is essentially the square 
root of the ratio of the centrifugal force to an inertia force. This 
can be seen from the following transformation: 


_ JL Pr^r^^ 

p2w2 ■ 2r^ /w2\ 

The right-hand term is, in addition to the length ratio L/2rj., the ratio 

of the centrifugal force per unit fluid voliome at the passage entrance 
to the inertia force that world be necessary to slow down the fluid par- 
ticle from an initial velocity Wj. to the velocity zero over a length 
equal to the channel length L. 

Use of these substitutions in equations (5), (lO), and (l2) gives 
the following dimensionless equations: 


^,2 


1 - 


rM2^2x> 


djtj 

w 




r „ ^ 2rM2pT dn (f L , , dT , T da 

(cos p) — + d|+a di 


(14) 



2K,dim(«l - 4T 




^ da 
a d4 


1 ^ 
a A^ 


‘^K,dim 




(15) 


(16) 


vhere 


^K,dim - 



(17) 


CALCULATION EROCEDUEE 

The two equations (5) and (12) or (14 ) and (16) can be used to 
calculate two of the parameters that depend on x if the rest of the 
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parameters and an appropriate set of boundary conditions are prescribed. 

Two problems are most frequently encountered in applying this analysis 
to design calculations: (l) A passage of prescribed geometry is ..closed 

at one end and fed at the other end by gas of a prescribed state. The 
pressure along the outside of the porous passage wall of prescribed 
material (Cj^ and n known) and the local distribution of the flow ejec- 
tion through the porous wall axe prescribed as well. The local distribu- 
tion of the permeability of the channel wall has to be determined in such 
a way that the required flow rate through this wall (for instance, to ob- 
tain a desired wall temperature distribution) is obtained. (2) The passage 
geometry as well as the permeability of the porous portion of the channel 
walls of known material is prescribed. The channel is again closed at 
one end and fed at the other end with gas of a prescribed state. The 
outside pressure distribution is known. The distribution of the local 
flow rate through the porous wall has to be determined. 

The first problem is the easier one to calculate, because, with 
the local flow rate prescribed^ the mass flow within the channel may be 
obtained from integration of equation (6). Equation (s) or (l4) then 
can be used to calculate the internal pressure distribution along the 
passage, because all the other terms in either equation are known, and 
equation (lO) (or (l5)) can then be solved for the constants Cj<;; or 

dim.^ which determine the required local permeability of the porous 
wall. 

Case (2) mak:es it necessary to solve equations (5) and (l2) or (14) 
and (l6) simultaneously for the two unknowns pj^ or and m/A or 

|i, the internal pressure and the mass velocity, respectively. 

A simplification of the outlined calciilation procedure, investigated 
for assumed radial passages, considers only centrifugal and pressure 
forces in the balance of forces acting on the control area. Equation (l) 
reduces to 


The mamerical examples presented later in this report are obtained for 
an assijmed linear vaxiation in cooling-air temperature) as a consequence, 
T also varies linearly in the form 



and equation (14) reduces to 



(18) 


T = 1 + aC 
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With this expression for t inserted in equation (l8)^ direct integra- 
tion results in 


’^i,s 


ii! (— 

2 I L 


In (l 




(19) 


The dimensionless pressure ^ in the passage is now independent of 

the coolant flow. The conditions under which the simplified procedure 
can be used are discussed later. 




NUMERICAL EXAMPLES 

Sample calculations for each of the two problems previously discussed 
were made for cooling-air passages I and II of the turbine blade shown 
in figure 2. The parameters needed in the calculations are indicated in 
the figure legend. They were chosen to correspond to conditions prevail- 
ing in a test engine available at the NACA Lewis laboratory for the test- 
ing of turbine blades. Passages I and II were assumed to be radial pas- 
ages (p = O). The cross-sectional area of passage II is assumed to 
decrease linearly along the blade span to a value at the tip which is 
half that at tne root. That of passage I is constant for all calculations 
except one^ the case where the effect of area variation on flow conditions 
was investigated. For this case^ two arbitrary passage area variations 
to be discussed later were considered^ in addition to the case of constant 
area. Poroloy wire cloth was selected as the porous material for the 
turbine blade skin^ constituting the porous side of each passage. An 
empirical equation supplied by the manufacturers of Poroloy (ref. 5) and 
valid for flow rates between 0.0001 and 0.1 pound per second per square 
inch gives the value n in equation (lO) as 5/8. The expre^^sion 

will be used later to determine the permeability per unit thickness K/t 
of the porous side of passages I and II (see fig. 2). 

The external pressure variations along passages I and II are also 
shown in figure 2. These spanwise pressure distributions were obtained 
from calculated velocity and pressure distributions around the blade 
periphery at three spanwise locations. The effective gas temperature 
profile along the blade span shown in figure 2 Vas determined experimen- 
tally with uncooled blades of the same geometry in a test engine. The 
inlet cooling- air temperature was assumed as Tp = 640^ R. A linear in- 
crease of approximately 200^ F in cooling- air temperature through the 
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blade span was assumed arbitrarily. The friction factor f indicated 
in figure 2 was obtained from the relation valid for turbulent flow 
through a tube. It corresponds to a Reynolds number of about 2000 when 
this parameter is based on the hydraulic diameter of the passage. Ac- 
tually, the Reynolds number, and therefore the friction factor, varies 
along the passage. This variation was neglected in the present 
calculations . 

The following dimensionless quantities were obtained from the pa- 
rameters in figure 2; they are held fixed throughout the calciilations: 

f = 0.0475 
T = 1 + 0.328 K 
r^/L = 2.569 

K = 0.2916 


Permeability Requirements for Prescribed Flow Ejection 

The determination of the permeability of the porous wall necessary 
to fulfill prescribed flow ejection distributions for the two passages 
of the blade shown in figure 2 will be considered first. The results 
of these calculations are presented in figures 3 to 6. The two pass^es 
were chosen so that laminar boundary- layer flow is expected along the 
outside blade surface adjacent to one passage (I) and turbulent flow 
along the outside blade surface adjacent to the other passage (ll) when 
transition is estimated from the pressiire distribution around the blade 
periphery. A blade wall temperature of 600^ F was prescribed. The 
spanwise cooling-air-flow rate pV can be calculated from the data pre- 
sented and from the condition that the blade temperature is constant over 
the surface (ref. 4). The cooling-air-flow distribution and a knowledge 
of b^ the circumferential length of the porous section Of the passage 
at each spanwise position, permit the integration of equation (6) and 
hence the determination of the distribution of the mass flow m through- 
out the passage. After conversion to the dimensionless mass-flow pa- 
rameter [I, equation (l4) can be solved for the dimensionless internal 
pressure distribution For the integration of equation (14), the 

cooling-air pressure at the entrance of the passage must be known. For 
this calculation, the value of the dimensionless external pressure at 
the root was assumed to be 0.95. From the internal and external 

dimensionless pressure ratios 7i^ and respectively, the dimension- 
less permeability parameter is obtained from equation (15). 

This parameter is valid for any porous material for which the mass flow 
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Is described by eqiaation (12). For Poroloy the desired permeability per 
unit thickness is determined from the expression 


K 

t 



1 2 -- 

, n / pv \ n 1 T 

•K.dim] p2.^To 


resxilting from use of equations (17) and (ll). 

Figure 3 shows ^ in dimensionless form^ the results of this calcula- 
tion^ namely^ the internal mass-flow distribution |i (which is nearly 
linear) and the internal pressure distribution as functions of the 

dimensionless distance C = x/L for passage whose porous wall is 
exposed to laminar external flow. The prescribed external pressure 
distribution is also shown once more for comparison. Figure 4 shows 

the dimensionless permeability parameter which is valid for any 

porous material for which equation (lO) holds ^ and the permeability dis- 
.tribution K/t required to satisfy the prescribed cooling-air ejection 
rates through the porous Poroloy wadi. A rapid reduction in blade per- 
meability or a corresponding increase in wall thickness is required near 
the lower quarter- span of the passage^ a more gradual permeability re- 
duction or thickness increase to about the three-quarter- span location^ 
and a nearly constant permeability or thickness in the region near the 
blade tip. Such a permeability or thickness variation will be extremely 
difficult to fabricate; and, in all probability, some sort of compromise 
between prescribed wall temperature (and corresponding prescribed 
coolant-flow ejection pV) and a more readily obtainable permeability or 
thickness variation will be necessary. 

A variation of the parameter K/t as shown in figure 4 is also often 
unfavorable with regard to stress, because it may mean heavier material 
for the blade shell near the tip of the blade. It is easily understand- 
able that the variation of the parameter K/t along the blade length be- 
comes smaller when the pressure at the passage inlet is increased. This 
is shown in more detail in some later figures. The result of the simpli- 
fied calculation procedxire (eq. (19 )) is shown in figure 3 as the dashed 

line jtn- o- It can be observed that the difference between the result 
X, s 

of the original calculation (the solid line itj^) and the result of the 
simplified procedure (the dashed line ^i^s) very small under the as- 
sumed parameters. 

Passage II differs from passage I by the fact that the gas flow 
around the blade is expected to produce a turbulent boundary layer adja- 
cent to passage II. Consequently, the coolant-flow ejection pV 
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required to maintain a blade wall temperature of 600^ F is considerably 
larger than for passage I. This distribution along the span of passage 
II was calculated by the method discussed in reference 4 and is shown in 
figure 5. The cross-sectional area of passage II was assumed to decrease 
linearly from blade root to tip; its value at the tip was half its value 
at the root. In algebraic form^ the dimensionless area of passage II 
may be expressed as a = 1 - 0.5 The calculation proceeded as in the 
previously described case. The results are plotted in figures 5 and 6. 
Dimensionless distributions of external pi^essure internal pressure 

and mass velocity p are shown as functions of the dimensionless 

distance C = x/l in figure 5. The variation of the cross-sectional 
area of the passage is reflected in the distribution of mass velocity 
m/A. The dimensionless permeability parameter C^r and the distri- 

bution of K/t are presented in figure 6. The variation of these pa- 
rameters over the blade height is qualitatively the same as for passage 
I. The required permeabilities^ however^ are considerably leirger for 
passage II in order to accommodate the lajrger coolant-flow rates. 

The dimensionless internal pressure distribution „ calciilated 

by the simplified equation is inserted in figure 5 as the dashed line. 

The difference between the solid line pressure distribution obtained 

by use of equation (l4) and the dashed line in figure 5 is larger than 
it was for passage I (fig. 3). Since the calciilation is considerably 
simplified by use of equation (18) instead of equation (14)^ it is im- 
portant to know when this simplification is permissible. This will be 
discussed in detail later. 


Flow Ejection Distribution for Prescribed Permeability 

The determination of the flow ejection distribution for a prescribed 
porous- wall permeability requires the simultaneous solution of equations 
(14) and (16). The two boundary conditions that must be fulfilled apply 
at opposite ends of the passage; that is^ the dimensionless internal 
pressure at the passage entrance ^ must equal unlty^ and the dimen- 
sionless mass velocity \i must reduce to zero at the passage tip. More- 
over^ must exceed (or in the limiting case^ equal) the dimensionless 

external pressure everywhere aJLong the passage to ensure flow through 

the wall in the proper direction. For a prescribed distribution of the 
outside pressure a unique p distribution that reduces to zero at 

the passage tip can be found. However^ this distribution can only be 
found by a trial- and- err or procedure. It is necessary to assume values 
of Pj^ p and (m/A)p at the passage entrance, to solve the equations 

numerically, and to check the value of p at the passage tip. In gen- 
eral, several trials will be required before the value p = 0 at the 
tip is obtained. 
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Calculations were made for the determination of the flow ejection 
distributions in passage I for three assumed values of the internal pres- 
sure at the passage entrance pj_ ^ (and hence for three distributions of 

^e “ Pe/Pi r) ^ passage with a porous wall whose permeability is 

Ck = 10"^ ftV2/( sec )(ibV4). Uniform porous-wall thickness and constant 

passage area (a = l) were assumed. The external pressure distributions 
jtg, the internal pressure distributions and the internal flow dis- 

tributions |i axe shown in figure 7. Also shown, as a dashed curve, is 
the internal pressure distribution g, which includes rotational ef- 
fects only and was calciilated by use of equation (19). Only slight varia- 
tions in the internal pressure distribution Jtj[ resxilt from variation 

of the press\ire with which the cooling air enters the blade root. More- 
over, the Internal pressure distribution obtained by use of equa- 
tion (18) is only slightly different from the other distributions 

obtained by use of eq-uation (14). For this blade passage, therefore, the 
simplified procedure wo-uld give practicably the same resiolts as the more 
complicated solution Involving equation (l4). The saving in time is 
especially large in this case, since no trial- and- error procedure has to 
be employed when equation (18) is used. 

The distribution of the air ejection rate pV is shown in figvire 8. 
As expected, the pV values increase as «e,r 'ieci’eases. Also, the 

ratio of the pV value at any location to the value at the passage en- 
trance changes less along the blade length with increasing inlet pressure. 
The pV distribution that results in a constant blade shell temperature 
is nearly constant spanwise, as can be seen from figure 3. A constant 
value of K/t for the porous wall will therefore lead to a blade tem- 
perature that decreases from root to tip. To make this decrease small, 
the cooling-air inlet pressure at the blade root sho\ild be made as large 
as possible. 

Figure 8 can be used to estimate how much more cooling air is re- 
quired for passsige I when the porous wall has constant permeability than 
when the permeability varies spanwise so as to give the constant or nearly 
constant coolant-flow rate pV necessary to maintain a constant spanwise 
wall temperatiire. For purpose of comparison, it may be assianed that the 
wall permeability and inlet pressure at the passage entrance are adjusted 
so that the coolant flow pV (and hence also the wall temperature) at the 
base is the same for both cases. It is then seen that, for the cases of 
a constant wall permeability with =0.8 and 0.6, the corresponding 

excess use of cooling air over the ideal case of a constant spanwise pV 
is approximately 77 percent and 45 percent, respectively. For the limit- 
ing case where Jtg p = 1.0, the coolant flow pV is zero at the blade 

root, and the condition of a constant spanwise coolant flow pV reduces 
to the case of an impermeable blade wall. 
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A series of calculations was made for the same passage^ the same 
wall permeability^ the same external pressure distribution corresponding 
to ^e^r ” 0*8^ but for passages whose flow areas increase (a =1 + 0.5 

are constant (a = l)^ and decrease (a = 1 - 0.5^ ) in flow direction, 
respectively. Figures 9 and 10 present the res-ults of these calculations. 
Figure 9 shows very little variation in the various distributions 

(including spite of rather marked variations in the mass ve- 

locity distributions \i. Accordingly, the pV distributions indicate 
(fig. 10) that passage flow-area variatibn has an almost negligible in- 
fluence on the flow ejection rates. 

From the fact that the curves in figure 9 are all in reasonably 

close agreement with the g curve resulting from the simplified cal- 
culation (eq. (19)), it may be deduced that changes in the variation of 
cross-sectional area do not cause appreciable differences between the 
pressure distributions obtained by the two calculation procedures. It 
can be concluded, therefore, that the rate of coolant flow, or Mach num- 
ber, is the determining parameter in this regard. 

Examination of equations (l4) and (18) reveals that the two solutions 
are identical for the case of zero coolant flow (M = 0)j and, hence, for 
small Mach numbers M good agreement .between the solutions for and 

^i,s (14) and (l8), respectively) is expected. That this is the 

case for the present calculations can be seen from figures 3, 5, 7,^ and 
9. In figure 3, the Mach number is very small (M = 0.059), and the 
solutions for and g are practically identical. In figure 7, 

the solutions for are for Mach numbers M of 0.135, 0.185, and 

0.228j the corresponding maximum percentage differences between and 

JCf “ g 

jti g are ^xl00% = 0.9, 1.8, and 3.1 percent, respectively. In 

^ ^i,s 

figure 9, where all solutions are for an essentially constant Mach number 
of 0.185, the maximum percentage difference between and g is 

only 2.1 percent. The solution for presented in figure 5, however, 

represents a comparatively high coolant-flow rate with a Mach number of 
0.350, and the maximum percentage difference between and is 

8.6 percent, an appreciable amount. From these results, it appears that 
the use of the simplified solution for the internal, pressure ^ given 

±, b 

by equation (l9) resiilts in an error of order not more than approximately 
2 percent from the exact solution, as long as the coolant-flow Mach num- 
ber does not greatly exceed 0.2. 
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SUMMARY OF RESULTS 

A one-dimensional method that permits the calculation of local per- 
meability for prescribed coolant flow (or the solution of the inverse 
problem) for a rotating passage with varying cross section and partially 
porous walls was developed. A knowledge of the temperature distribution 
of the cooling air inside the passage, the outside pressiire along the 
passage, and the passage geometry is required. Solutions to the first 
problem were obtained by solving independently a pair of first-order dif- 
ferential equations. For the inverse problem, a simultaneous solution 
of a pair of equations is required) a trial- and- error procedure has to 
be applied in this case. 

Solutions to both problems were obtained for several blind radial 
passages in a transpiration- cooled turbine rotor blade. Tlie results of 
this investigation are as follows: 

1. For prescribed coolant-flow ejection rates necessary to maintain 
a constant porous-wall temperature for both a constant-area passage ex- 
posed to external laminar-flow conditions and a variable-area passage 
exposed to external t\arb\ilent-flow conditions, considerable reduction in 
permeability or increase in thickness is required from blade root to 
blade tip. In either case, a rapid reduction in permeability or increase 
in thickness is required near the root of the blade) the reduction tends 
to level off considerably beyond this point. 

2. Permeability or thickness variations like those necessary to 
maintain a constant wall temperature throughout appear extremely diffi- 
cvilt to obtain with present-day materials, and some compromises between 
prescribed wall temperature and a more readily obtainable permeability 

or thickness distribution will undoubtedly be necessary. These variations 
will probably require some increase in shell thickness near the tip of 
the blade and, as a consequence, they might- be unfavorable from stress 
considerations . 

3. For prescribed locally constant permeability, the mass-flow ejec- 
tion rate increases from blade root to blade tip. As the passage inlet 
pressure is increased, the relative increase in flow ejection rate along 
the blade height becomes less. Consequently, a large inlet pressvire will 
result in a more nearly uniform flow ejection rate. 

4. For prescribed locally constant permeability, variations in pas- 
sage cross-sectional area had practically no influence on flow ejection 
rates . 

5. The internal press-ure distribution of cooling air in a rotating 
passage can be approximated by balancing the centrifugal force due to 
rotation with the pressure forces only. This approximation can be 
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applied to the solution of both types of problems discussed above. For 
many cases of practical importance^ the simplified differential eqiiation 
can be Integrated directly to obtain the pressure distribution^ thus 
eliminating the numerical process necessary for the complete solution. 

The complete solution and the approximate solution converge for zero 
coolant flow. The res\ilts presented in this report indicate that the 
complete solution for the internal pressure is well approximated by the 
simplified solution, for cooling-air Mach numbers at the passage entrance 
^ that do not greatly exceed 0.2. The maximiim difference between the com- 

^ plete and approximate solutions is then of the order of 2 percent. 

CO 
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Figure 1. - Transpiration- cooled turbine blade and detailed element of 
control area used in derivation of basic equations. 
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Figure 2. - Spanwlse distributions of effective 
gas temperature and external pressure. Angular 
velocity, o), 1205 radians per second; length, L, 
0.3 foot; radius, r^, 0.7708 foot; coolant tem- 
perature, T, 640°+700x (^); friction factor, 
f, 0.0475. 
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Figure 3. - Prescribed external pressure and coolant-flow ejec- 
tion distributions and resulting internal pressure and coolant 
mass-velocity distributions for passage I. 
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Figure 4. - Permeability distribution required to satisfy pre- 
scribed coolant-flow ejection distribution for passage I. 
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Figure 5. - Prescribed external pressure and coolant-flow ejec- 
tion distributions and resulting internal pressure and mass- 
velocity distributions for passage II. 
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Figure 6. - Permeability distribution required to satisfy pre- 
scribed coolant-flov ejection distribution for passage II. 
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Figure 7. - Effect of 
prescribed constant 


variation of passage Inlet pressure on Internal pressure and mass-velocity distributions for passage I with 
wall permeability Cj^ » 10"^ ft^/^/(sec 
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Figure 8. - Coolant-flow ejection distributions for passage I 
for various Inlet pressures and prescribed permeability 

Ck = 10"^ sec) (ibV4) . 
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(a) External pressure. 


(b) Internal pressure. 


(c) Mass velocity. 


Figure 9. - Effect of passage area variation on internal pressure and mass-velocity distributions for prescribed external pressure 
distribution and wall permeability Cj^ = lO”^ ft^/^/(sec) (lbV4=) . External pressure ratio, ~ Pe,r/Pi,r “ 0.8. 




NACA TN 3408 


NACA TN 3408 


29 



Figure 10. - Coolant-flow ejection distributions for passage 
with prescribed area variations, external pressure distribu- 
tion, and wall permeability = 10”^ ft^/^/(sec )(lW^). 
External pressure ratio, «e,r = Pe,r/Pi,r 0 , 8 . 
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